Field theory of bicritical and tetracritical points. IV. Critical dynamics including 

reversible terms. 
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This article concludes a series of papers (R. Folk, Yu. Holovatch, and G. Moser, Phys. Rev. E 
78, 041124 (2008); 78, 041125 (2008); 79, 031109 (2009)) where the tools of the field theoretical 
renormalization group were employed to explain and quantitatively describe different types of static 
and dynamic behavior in the vicinity of multicritical points. Here, we give the complete two loop 
calculation and analysis of the dynamic renormalization-group flow equations at the multicritical 
point in anisotropic antiferromagnets in an external magnetic field. We find that the time scales 
of the order parameters characterizing the parallel and perpendicular ordering with respect to the 
external field scale in the same way. This holds independent whether the Heisenberg fixed point or 
the biconical fixed point in statics is the stable one. The non-asymptotic analysis of the dynamic 
flow equations shows that due to cancelation effects the critical behavior is described - in distances 
from the critical point accessible to experiments - by the critical behavior qualitatively found in one 
loop order. Although one may conclude from the effective dynamic exponents (taking almost their 
one loop values) that weak scaling for the order parameter components is valid, the flow of the time 
scale ratios is quite different and they do not reach their asymptotic values. 



I. INTRODUCTION 



Three component antiferromagnets in three spatial di- 
mensions in an external magnetic field in z direction con- 
tain in their phase diagram two second order transition 
lines: (i) between the paramagnetic and the spin flop 
phase and (ii) between the antiferromagnetic and para- 
magnetic phase. The point where these two lines meet is 
a multicritical point which turned out to be either bicrit- 
ical or tetracritical. Within the renormalization group 
(RG) theory the stability and attraction region of the 
static fixed point (FP) of the RG transformation deter- 
mines, which kind of multicritical behavior is realized. 
For the bicritical point it is the Heisenberg FP, for the 
tetracritical point it is the biconical one. The stabilty of 
a FP depends on the system's global features as the space 
and order parameter (OP) dimensions d and n. In d = 3, 
the case considered here, the biconical FP is stable apart 
from a restricted attraction region of the Heisenberg FP 
. The static phase transition on each of the phase transi- 
tion lines belongs for (i) to an isotropic Heisenberg model 
with n± = 2 and for (ii) to Heisenberg model with = 1 
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Concerning the dynamical universality classes the 
transition (i) belongs to the class described by model F 
and (ii) belongs to the model C class (for the notation 
see Q). At the multicritical point the critical behavior 
is described by a new universality class both in statics 
and dynamics. The interesting feature of these systems 
is that all the different OPs characterizing the ordered 
phase are physically accessible. This is most important 
for the dynamical behavior since the only other example 
belonging to model F is the superfluid transition in 4 He 
where the OP (the complex macroscopic wave function 
of the condensate [i| ) is experimentally not accessible Q ■ 
Here the OPs are the components of the staggered mag- 
netization. Their correlations (static and dynamical) can 
be measured by neutron scattering. 

A complete description of the critical dynamics near 
the multicritical point mentioned above has to take into 
account the slow dynamical densities which are the OPs 
and the conserved densities present in the system. Due 
to the external magnetic field the only conserved density 
which has to be taken into account is the magnetization 
in direction of the external field. A derivation of the 
dynamical equations follows along the usual steps cal- 
culating the reversible terms from the non-zero Poisson 
brackets, introducing irreversible terms present also in 
the hydrodynamic limit, dropping irrelevant terms and 
taking into account terms arising in the renormalization 
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procedure (see e.g. the review Q). Such a dynamical 
model has already been considered in by RG the- 

ory and it was argued that due to nonanalytic terms in 
e = 4 — d a FP in two loop order qualitative different 
from the one loop FP is found. The result of the one 
loop calculations is that the time scales of the parallel 
and perpendicular components of the staggered magne- 
tization scale differently whereas calculated in two loop 
order they scale similar although the FP value of the 
timescale ratio of the two components cannot be found 
by e expansion and might be very small in d = 3 namely 
of 0(1O~ 86 ). It was argued that the terms leading to 
the singular behavior in e do not contribute to the FP 
value of the mode coupling. The calculations of the RG- 
functions in where not complete in two loop order 
(they took into account only the terms which lead to the 
nonanalytic behavior in e). At that time also the Heisen- 
berg FP (named TL) was considered to be the stable static 
one, whereas it turned out in two loop order (resummed) 
that it is the biconical FP (named B) [2|. 
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TABLE I: Dynamical FP values (zeros of the corresponding 
dynamical /3-functions) at d = 3 of different models for the 
time scale ratios wt , w'2 , v* and the mode coupling constant 
fj_ . The 2nd and 3rd lines quote results of this paper found in 
the biconical static FP for the tetracritical behavior of the dy- 
namical model that takes into account reversible terms. They 
are compared with the two loop results found in the model 
C multicritical point [ll| as well as in the critical points of 
model C for the one component OP [9] and of model F for 
the two component order parameter [10j . 

A summary of the results obtained so far for the FPs 
characterizing dynamical behavior is given in table HI Ne- 
glecting the reversible terms one is left with a purely 
relaxational dynamics. Then the asymptotic dynamical 
critical behavior is characterized by the FP values of the 
independent time scale ratios of the system. These are 
the following time scale ratios: (1) the ratio w\\ between 
the relaxation rate of the staggered magnetization paral- 
lel to the external field and the diffusive transport coef- 
ficient of the magnetization parallel to the external field; 
(2) the ratio w' ± between the real part of the relaxation 
rate of the staggered magnetization perpendicular to the 
external field and the diffusive transport coefficient of the 
magnetization parallel to the external field. In addition 
we introduce the ratio v' between the two components 
of the real relaxation rates of the two OPs in order to 
compare their dynamic scaling behavior. A non-zero fi- 
nite value of the time scale ratio means that the two 
involved densities scale with the same exponent. If all 



time scale ratios are non-zero and finite, one speaks of 
strong dynamic scaling, otherwise of weak dynamic scal- 
ing. Especially of interest is the behavior of the scaling 
of the two components of the OP indicated by the FP 
value of v' . In the third line of table U the two loop order 
result shows weak dynamic scaling between the OPs and 
the conserved density but strong scaling between the OP 
components. However since the FP value time scale ratio 
v' is almost zero, the critical behavior is dominated by 
non-asymptotic effects. For comparison the FP values 
for the case of model C for the one component OP @ 
and for model F for the two component order parameter 
A (j are included. They are the limiting cases when the 
two OPs characterizing the multicritical behavior would 
decouple in statics and dynamics. 

In the first line of table U the results for the multi- 
critical dynamical FP Be values taking into account the 
static coupling of the OP to the conserved density are 
displayed (see [HI). All time scale ratios are non-zero 
and finite but since w'^_ is very large (v 1 * almost zero) 
the observable behavior in the vicinity of the multicriti- 
cal point is predicted to be dominated by non-asymptotic 
effects and strong scaling is not observable [11(. In the 
second line the results of a one loop RG calculation with 
reversible terms for the biconical FP are given. The FP 
value of the mode coupling parameter f± is finite but 
since wt, = the critical dynamics is characterized by 
weak dynamic scaling and the two components of the 
OP scale different. A similar result for the Heisenberg 
FP was found in In the third line the results found 
in this paper are shown, indicating weak scaling between 
the conserved density and the components of the OP, 
but strong scaling between the parallel and perpendicu- 
lar components of the OP. Since the FP value of the time 
scale ratio between the component v'* is almost zero but 
definitively different from zero it is expected that non- 
asymptotic behavior is dominant. 

This article concludes a series of papers @, [n]> [l2| 
(henceforth cited as papers I, II, and III) where the tools 
of the field theoretical RG were employed to explain and 
quantitatively describe different types of static and dy- 
namic behavior in the vicinity of multicritical points. A 
short account of the results presented here was given in 
[r|. The statics and dynamics were treated in Refs. Q 
and [HI, respectively. First, purely relaxational dy- 
namics was considered (paper II) and later, in paper III, 
these results served as a basis to consider how an account 
of magnetization conservation modifies dynamical behav- 
ior. The goal of the current study is more ambitious: we 
will analyze a complete set of dynamical equations of 
motion taking into account reversible terms [l4|, [H| and 
give a comprehensive description of dynamical behavior 
in the vicinity of multicritical points in two loop order. 
The paper is organized as follows: In section [IT] the dy- 
namic model is defined followed by a the definitions of 
the dynamical functions considered in section Mil The 
renormalization and corresponding RG- functions are pre- 
sented in section ITVl and fVl respectively. The two loop re- 
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suits of our calculations for these dynamic RG-functions 
are given in section fVTl The one loop approximation for 
the dynamic is discussed in section fVIII In the next sec- 
tion IVIIII we consider the asymptotic properties of the 
two loop RG-functions leading to the general asymptotic 
results in section IIXI We then present the results ex- 
pected in the asymptotic subspace, section [X] The non- 
asymptotic behavior, obtained by looking at the region 
further away from the multicritical point, is shown in sec- 
tion [XI] a conclusion IXIII ends the paper. In Appendices 
calculational details for some intermediate steps of the 
RG calculation are presented. 



II. MODEL EQUATIONS OF THE 
ANTIFERROMAGNET IN AN EXTERNAL 
FIELD 

The non-conserved OP 4>o of an isotropic antiferromag- 
net is given by the three dimensional vector 




(1) 



of the staggered magnetization, which is the difference 
of two sublatticc magnetizations. An external magnetic 
field applied to the ferromagnet induces an anisotropy to 
the system. The OP splits into two OPs, (j>±o perpen- 
dicular to the field, and <fi\\ parallel to the external field. 
Assuming the z-axis in direction of the external magnetic 
field, the two OPs are 




(2) 



In addition the z-component of the magnetization too has 
to be taken into account for the dynamics and therefore 
has to be included in statics although there it could be 
integrated out and does not change the asymptotic static 
critical behavior. Thus the static critical behavior of the 
system is described by the functional 
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H= I d xi -f±(f>±o ■ 4>ia , 
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with familiar notations for bare couplings {u,j}, masses 
{r} and field h The critical dynamics of relaxing 



OPs coupled to a diffusing secondary density is governed 
by the following equations of motion [6|: 
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(4) 



(5) 
(6) 



with the Levi-Civita symbol e a ^ z . Here a, (i = x,y and 
the sum over repeated indices is implied. 

The dynamical equations describe the dynamics of an 
antiferromagnet with the usual Lamor precession terms 
for the alternating magnetization and relaxational terms. 
Due to the static coupling to the conserved magnetiza- 
tion additional Lamor terms arise together with a diffu- 
sive term for the magnetization. Renormalization con- 
siderations lead on one hand to a neglection of several 
Lamor terms and on the other hand create an additional 
reversible term (the second term on the right hand side 
of Eq. [5])) not present in the usual dynamics of antifer- 
romagnets [16j . 

Combining the kinetic coefficients of the OP to a com- 
plex quantity, T±_ = r j_ + IT" , the imaginary part con- 
stitutes a precession term created by the renormalization 
procedure even if it is absent in the background. The ki- 
netic coefficient A and the mode coupling g are real. The 
stochastic forces 6 ( j >± , 6^,„ and 9 m fulfill Einstein relations 

(6$ ± (x,t) 9^(x',t'))=2t' ± S(x-x')S(t-t')6^ , (7) 
(6 H {x,t) 9 H (x',t'))=2f ll S(x-x')S(t-t r ) , (8) 
(6 m (x,t) 6 rn {x' ,t')) = -2\V 2 8{x -x' )5{t-t') . (9) 

In view of dynamical calculations it is more convenient 
to deal with a scalar complex order parameter ipo = ip' + 
iipQ instead of the real two-dimensional OP <fij_o in ((2). 
Thus we may introduce 



(10) 



as OP of the perpendicular components. The superscript 
+ denotes complex conjugated quantities also in the fol- 
lowing equations. In addition to the two OPs the z- 
component of the magnetization, which is the sum of the 
two sublattice magnetizations, has to be considered as 
conserved secondary density Too- 

Expressed in terms of the above densities the dynamic 



4 



equations take the form 



extended static functional (|16|) by 
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(11) 
(12) 
(13) 
(14) 



Due to the fact that the stochastic forces 9% in (j4]) are 
^-correlated and fulfil the Einstein relations, similar prop- 
erties hold also for the stochastic forces 9^: 

(9^(x,t) 9+(x',t'))=4f' ± 6(x-x')6(t-t'). (15) 

The critical behavior of the thermodynamic derivatives 
follows from the extended static functional (the func- 
tional ((3]) written in the variables intoduced in (QUI)) 



n = n {0) + n [mt) (i6) 

with a Gaussian part 

H<® = J d d xU^+^ + I(V^ +)(V^o) 

+^||^o + |(V^|| ) 2 + \rnl - Am | , (17) 
and an interaction part 
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The above static functional may be reduced to the 
Ginzburg-Landau- Wilson (GLW) functional with com- 
plex OP by considering an appropriate Boltzmann dis- 
tribution and integrating out the secondary density. One 
obtains 

H GLW = J d d x\U^+^ + I(W + )(V^ ) 



2 """ u 2* 



(19) 



The parameters {r} = f±,rn and {u} = iix,ui\,u x in 
([in]) are related to the corresponding parameters of the 



fx = f± + 7±/i , 
f|| = f|| +7||ft , 
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(20) 
(21) 
(22) 



The property that the static critical behavior does not 
depend on the secondary densities, which can be inte- 
grated out in (TTB1) . leads to relations between the corre- 
lation functions of the secondary densities and the OP 
correlation functions. These relations and their deriva- 
tions have been extensively discussed in paper III with 
real OP functions </>xo and 4>\\o- Because the derivation 
of the relations is independent of the type of OP (real 
or complex) , all of the relations remain valid and can be 
taken over from paper III. Therefore we will not repeat 
them here. 



III. DYNAMIC CORRELATION AND VERTEX 
FUNCTIONS 

The Fourier transformed dynamic correlation functions 
of the two OPs are usually introduced as 



CW+ ({£}> k ,u) = d d x / dte 



-ikx+iojt 



(Vb(a:,t)l#(0,0)) e 
(23) 



CV*n*i|({£;}> fc > w ) = d d x dte 



^ + ^(0||o(^,O0||o(O,O)) c 
(24) 

All functions depend on the two correlation lengths £x 
and which is indicated by {£} in a short notation. 
{AB) C = (AB) - (A)(B) denotes the cumulant. The av- 
erages are calculated with a propability density includ- 
ing a dynamic functional, which can be constituted from 
the dynamic eq uations (fTTj) - (fT4"|) . In the considered 
approach of [17j | for every density auxiliary densities are 
introduced accordingly. They are denoted as ipQ , tpO: 4>\\o 
and rfiQ. The dynamic correlation functions of the order 
parameters are connected to dynamic vertex functions 
via 



CW+({£}' fc > w ) 



r w+ (ft},*,") 



f^ + (tf},fc, w ) 



(25) 



(26) 



where the two-point vertex functions appearing on the 
right hand side in the above expression have to be calcu- 
lated within perturbation expansion. They are obtained 
by collecting all 1-particle irreducible Feynman graphs 
with corresponding external legs. A closer examination 



5 



of the loop expansion reveals [18| that the dynamic re- 
sponse vertex functions and i have the general 
structure 

+ f^ + (U},*)f^ + ({C},*,w),(27) 



+ f 0||0|| ({e},fc)f||, (28) 

where f^^+ ({£}, fc) and f ^..^ ({£,}, k) are the well known 
static two point vertex functions of the bicritical GLW- 
model with a complex OP. We want to remark that the 
static vertex functions in (|27|) and (|28| are related by 



and 
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(29) 



(30) 



to the static vertex functions introduced in papers I— III 
for the model with real OPs. Thus the correlation lengths 
and £|| are now determined by 



d In t^+(k, {f}, {u}) 
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(31) 



(32) 



are purely dynamic functions. 



The explicit expressions of these functions are given in 
appendix S Eqs. (|A"Tj) - ([AH . They determine also 
the dynamic vertex functions and ^ in (|25[) 

and (|26p . A proper rearrangement of the perturbative 
contributions shows that the relations 



^+({£h fc > w ) 



-23? 



ipip 



(33) 



f ^„ ({£}, M) = -2T,|K [ft^ ({£}, *, w)J (34) 

hold. is the real part of the expression in the brack- 
ets. 

Analogous to (f2"3"]l and the Fourier transformed 
dynamic correlation function of the secondary density is 
introduced as 

C mm ({S},k,u) = fd d xfdte- ikx+iut (m (x,t)m (0,0)) c 

(35) 



The connection to the dynamic vertex functions is anal- 
ogous to the case of the OP Eqs. ([231 and 



r< (ft:~\ 1 \ r>?im({£}, k, io) . , 

C mm ({£}, k,LJ) = (3 6 ) 



fmm({^}) fc > w ) 



The dynamic response vertex function of the secondary 
density has the general structure 

r m m(U},fc,w) = -iujrL m m({£},k,Uj) 

+ fm m ({C},*)f^({C},*,«) (37) 

where r m m({£},fc) is the static two point vertex func- 
tion calculated with the extended static functional (fT6| 
which already has been introduced in paper III. A re- 
lation corresponding to ()33l) holds also for the dynamic 
vertex function of the secondary density. We have 



rr«m({£},&,w) 



-25ft 



(38) 



IV. RENORMALIZATION 

A. Static renormalization 

The renormalization of the GLW-functional (fT9|) has 
been extensively discussed in paper I. The only difference 
in the present paper is that we now have to renormalize 
the complex OP V'o instead of the real vector OP <j>±o- 
We introduce the following renormalization factor 



(39) 



where Z$ is a real quantity and identical to Z^^ in paper 
I taken at n± = 2 and nu = 1. This means 



Zii, — Za 



(40) 



The renormalization of the parameters r±, rii and the 
couplings u±, u\\, u x appearing in (|19p is given in paper 
I (see Eqs. (16), (17) and (5)-(7)). In all relations one 
has to replace Z^ ± by Z^. All renormalization factors 
remain valid if one sets n± = 2 and nu = 1. This is also 
true for the Z-matrix Z^i introduced in Eq.(lO) of paper 
I and the additive renormalization A({it}) defined in Eq. 
(15) of paper I. 

The renormalization of the parameters in the extended 
static functional p7J|) has been presented in paper III. As 
in the case of the bicritical GLW-model all Z-factors and 
relations between them remain valid if Za x therein is 
replaced by and if one sets n± — 2 and nu = 1 in 
explicit expressions. 
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B. Dynamic renormalization 

Within dynamics auxiliary densities tpo, (f>»Q and m 
corresponding to the two OPs and the secondary density 
are introduced [l7|- Instead of renormalization condi- 
tions we use the minimal subtraction scheme [nj as in 
the preceding papers II and III. The auxiliary density 
of the complex OP is multiplicatively renormalizable by 
introducing complex Z-factors: 



where Z^' contains only the poles of the k 2 derivative of 

^mm taken at zero frequency and wave vector modulus. 

The mode coupling coefficient needs no independent 
renormalization, so we simply have 



g = n e/2 Z m gA 



-1/2 



(48) 



The geometric factor Ad [20| already used in the static 
renormalization has been given in paper I Eq. (8). 



-7 1 / 2 7H 



(41) 



The complex renormalization factors in (|4ip fulfill the 
relation = Z~\ For the auxiliary densities of the 

single-component real OP and the secondary density the 
corresponding renormalization factors are introduced: 



'llo 



= zV\ 



m = Zmih . 



(42) 



where Z^ and Z^ are real. Within the minimal subtrac- 
tion scheme the Z-factors of the auxiliary densities of the 
non-conserved OPs Z^ + and Z^ are determined by the 

e-poles of the functions Cl^ + and introduced in 

(|27p and (f2"B")) . The corresponding function of the con- 
served secondary density Cl m f n in (|37|) does not contain 
new poles. Therefore one has 



Z~ 



(43) 



where Z m has been introduced in Eq. (30) in paper III. 
The kinetic coefficients renormalize as 



Zr„ T\i 



X = Zx\. 



(44) 



The renormalization of the complex kinetic coefficient T± 
in (J3U) leads to a complex Zr ± , while the other two renor- 
malization factors in (14"4"|) are real valued. Zr ± can be 
separated into 



V. RENORMALIZATION GROUP FUNCTIONS 

In order to obtain the temperature dependence of the 
model parameters, as well as the asymptotic dynamic 
exponents, the RG functions, which are usually denoted 
as £- and ^-functions have to be introduced. 

A. General definitions 



In order to simplify the general handling of the RG 
functions we will use the uniform definition 



dlriK 



(49) 



for all (^-functions in statics and dynamics. The deriva- 
tive is taken at fixed bare parameters, {aj} denotes the 
set of static and dynamic model parameters which in- 
clude the static couplings {u} and {7}, the mode cou- 
pling g, and all kinetic coefficients r_L,Ix,ri|,A. The 
(^-function £ ai is calculated from the renormalization fac- 
tor Z ai introduced in the previous section. Thus at may 
denote a model parameter from the set {aj}, a density 
4>x, <t>\\, m, or a composite operator M. The ap- 
proach of the model parameters cti(l) to their FP values 
in the vicinity of the multicritical point is determined by 
the flow equations with the flow parameter I 



7 i/2 7 -i/2 „(<Q 



(45) 



(<!) 



where Z-, 



namic function f ~ , 



r ^ contains the singular contributions of the dy- 
, appearing in (|2T[) . 



The dynamic equation ([131) for the OP <fru contains no 
mode coupling term. As a consequence only the kinetic 
coefficient F11 appears in the dynamic vertex function 



instead of a function T 
can write 



(d) 



Therefore Z 



id) 



7 _ 7V 2 7-I/2 

z r„ - Z, Z- 



1 and we 



(46) 



Using Eq. (j4"3")l the kinetic coefficient of the secondary den- 
sity renormalizes as 



(50) 



with /3-functions 



fl*({«i(0}) - Oi(0(- Ci + Ca 4 ({Oj(0})) (51) 

Cj is the naive dimension of the corresponding parameter 
cti obtained by power counting. For the static couplings 
itj_, Ux or itii the naive dimension Cj is equal to e, while for 
7_l or 7|| and the mode coupling g it is e/2 respectively. 
All kinetic coefficients, these are Tj_, T^, rii and A, are 
dimensionless quantities, which means Cj = 0. 

The flow equations ([50)1 have fixed points at the zeros 
of the /3-functions. The FP values of the model parame- 
ters {a*} are defined by the equations 



Zx 



(d) 



(47) 



^(K}) = o. 



(52) 
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The FP is stable if all eigenvalues of the matrix d(3 ai /dak 
are positive or possesses positive real parts. Starting at 
values {otj(lo)} at an initial flow parameter value Zo , the 
flow equations can be solved numerically. The asymp- 
totic critical values of the parameters are obtained in the 
limit I — > 0. If a stable FP is present the flow of the 
parameters has the property 



lim{a,-(/)} 



(53) 



A set of FP values determines all static and dynamic 
exponents. The static relations between ^-functions and 
critical exponents have been extensively discussed in pa- 
pers I and III. The dynamic exponents are related by 



Hi 



2 + £, 



2 + Cr n , z m = 2 + C x (54) 



to the dynamic £-functions (see Q). In ([54)) the short 
notation £*. = Ccti({ a j}) nas been introduced. In the 
non-asymptotic background region effective dynamic ex- 
ponents are defined as 



*i e//) (o 
4 e//) (0 



2 + C rl (K(0l) 

2 + Cr„(K(0) , 
2 + Ca(K(0) • 



(55) 
(56) 
(57) 



where the flow of the parameters is inserted into the £- 
functions instead of the FP values. The effective expo- 
nents depend on the flow parameter, or reduced tem- 
perature accordingly. Relation (|5"3"|) makes sure that the 
effective exponents turn into the asymptotic exponents 
in the critical limit, that is 



limz[ e//) (/) 
i^o fc 



with 



k =_L, 



now complex quantities. The ratios in Eqs. ([59l) and 
(|60p are of course not independent as shown by 
the equality in ((60]). The structure of the dynamic 
^-functions presented subsequently further implies 
the introduction of the complex ratio 



v± = 



Wj_ 



(ii) The mode coupling parameter 



F = 



(61) 



(62) 



The above ratio does not necessarily have a finite 
FP value. Thus it may be more appropriate to use 
the ratio 



fx = 



F 



(63) 



in several cases, especially in the discussion of the 
flow equations and the fixed points. 

The flow equations for the ratios defined above can be 
found from the £- and /3-functions introduced in the pre- 
vious subsection. From the definition of the parameters 
in ([59]). ([55]) and the renormalization flU]) and (gSJ) we 
obtain together with (|4"9"]) the flow equations 



dw± 
dw\\ 



= wx (Cr_L ~ Ca) 



l- 



dl 
df± 



w 



ll (Cpu-Ca) 



(58) = 
di 



w ± 



(64) 
(65) 
.(66) 



B. Time scale ratios and mode coupling parameters 



From (|64[) and (1651) follows immediately the flow equation 
for the ratio 



It is convenient to introduce ratios of the kinetic co- 
efficients or mode couplings, which may have finite FP 
values. The following ratios will be used in the subse- 
quent sections: 

(i) The time scale ratios between the order parameters 
and the secondary density 



1 1I = V ^ 



Cr 



(67) 



which has been defined in ([611]) . 

The remaining task is to calculate the explicit expres- 
sions of the dynamic functions Cr ± , Cry and Ca i n two 
loop order. 



U)J_ 



A 



El 
A 



(59) 



From this we may also define the ratio between ki- 
netic coefficients of the two order parameters 



_ r ll 



w 



(60) 



which already previously has been used in the bi- 
critical model A and model C. Note that in con- 
trast to the two models mentioned, w± and v are 



VI. DYNAMIC RG-FUNCTIONS IN TWO LOOP 
ORDER 

The perturbation expansion of the dynamic vertex 
functions and the structures therein are outlined in detail 
in appendix [S] The outcoming expressions for the dy- 
namic renormalization factors in two loop order are pre- 
sented in appendix [B] With these expressions at hand we 
are in the position to obtain explicit two loop expressions 
for the RG ^-functions as expressed in the following. 
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A. Dynamic (^-functions of the OPs 



Relation (|46| between the Z-factors implies the rela- 
tions between the corresponding ^-functions 



(68) 
(69) 



The static (^-functions Q$ — Q ( j )x has been presented 
Eqs.(20) in paper I. Inserting flST) and (|B2) into (|49|) 
and (|68l) we obtain the dynamic ^-function for the ki- 
netic coefficient of the perpendicular components as 



2 u\D\ 



w±(l + w±) 3wi(l + 
I D 2 



2w^(l + w 1 ) 2 

1 7H-P-L /ux 1 TjfPj 
21 + wi 1 3 2 1 + to 



+cS ) (M.«j..«) (70) 

where we have introduced the coupling 

D± = wx-yx - iF . (71) 
The functions Ax, B± and Xx arc defined as 

A ± = io±7_l(1 - xiLi) + iFx_xiLi - D±L (72) 

B x = ^i7i(l - 2x 1 L 1 ) + F 2 {2x-L x + L R ) 
+2w ± 7_lLF(1 + 2x^x 1 L 1 ) - 2L Dl 

° 2 ^ + ( l + 2^1nil±^)!l (73) 



1 + 10JL 



X i = 1 + In 



2w 
1 + v 



I + 2w_l 



1 + ^m^^ 
v J 2 + v 



with 



5+ +vx +x 2 + (x 2 + + 2v'i) 



Li 



— 3vj 



(74) 



(75) 



We have used the following definitions in the above ex- 
pressions: 



x± = I ± v± 



L = 2 In 



xx — 2 + v± , 



?— , Li = In 

1 + 1 + 2^- 



(76) 



(77) 



Cp ({u}, vx, v) is the £- function of the kinetic coefficient 
of the perpendicular components in the bicritical model 



A, but now with a complex kinetic coefficient r_i_. It 
reads in two loop order 



(W({u},vx,v) 



U ± (r r 1 

, u x ( T {y) 1 
+ 36l^ "2 



with 



L^.m (1 + V) 



2 2, 2(l + v) 
in ■ 



v(2 + v) v 2 + v 



(78) 



(79) 



The dynamic ^-function of the parallel component is ob- 
tained by inserting Eq.(21) of paper I and (|B3|) into (02) 
and ([59]), The result is 



J||7|| I t0||7|| 



I + w\\ 2 1 + w 



l|7|| I 1 — 3 In ■ 



™|| 7f 



ifi-ginl)-^ 

1 + W|| \2\ 3 / 1 + 10|| 

l + 2to,i (l + wii) 2 \ 
In ■ 



W||7|| 

1+Wll 





\ Tl ] 






Iw'J 



1 + ton 



1 + 2WH 



F 



r 2 



2w^(l+w||) 



(80) 



The functions T± and Ti are defined as 



Ti = Dj 



1 + In ■ 



1 + u 



iH (1+v) In 



(81) 



(1 + — In ■ 



1 



1 + v 



\v+-^-(l+v)j (v+vx(l+v)) 



«+^(l+«) 



(82) 



Cr^ U -L' v ) i s * ne C-f urLC ti orL of the kinetic coefficient 
of the parallel component in the bicritical model A. With 
a complex Tx it reads 



(83) 
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with 

(x) 



L 



In 



(l+„) (_L +V ) 



-(1+v) 



vv± In 



>ln 



(84) 



VII. CRITICAL BEHAVIOR IN ONE LOOP 
ORDER 

Although the one loop critical behavior of the consid- 
ered system has already been discussed in Q we want 
to summarize the results in order to compare it with the 
considerably differing results of the two loop calculation. 
In one loop order the (-functions (fTO]) . (|80|) and ([87]) re- 
duce to 



B. Dynamic (-functions of the secondary density 



With relation ([47]) we can separate the static contri- 
butions to the (-function (a. Thus we have 



(a = 2( m + C{ d) 



(85) 



By separating the static from the dynamic parts in the 
^-functions one can take advantage of the general struc- 
tures appearing in the purely dynamic (-function (^ as 
well as in the static (-function ( m . Inserting n± — 2 and 
nil =1 into relation (40) in paper III ( m can be written 
as 



1 



1 



2 "T± + 



(86) 



which is valid up to two loop order. From the diagram- 
matic structure of the dynamic perturbation theory fol- 
lows 



Ad) _ 



f( 1 + e ) 



(87) 



The real function Q contains all higher order contri- 
butions beginning with two loop order. Setting Q = 
in (|87[) reproduces the one loop expressions of this func- 
tion. The function Q in the dynamic (-function of the 
secondary density (|87|) has the structure 

Q = ~»[* 2 ] (88) 

from which immediately follows that it is a real quantity. 
Xi reads 



w\(l + w±) 



1 1+U>j 
— h In ; 

2 1 + 



(89) 



+D+(1 +w±)- (V{ m) 7 ± + ioxlF) W 



where we have introduced the definitions 



(m)^(m) 



(m) 



w 



(m) 



In 1 



WJ_ 



W 



(m) 



(90) 



(91) 



Note that X 2 coincides with the corresponding function 
in model F in 0,[2l|. 



Dl 



w±(l + w±) 



Ad) 
^A 



fl 



(92) 



Inserting (li?2l into the right hand sides of (|S4"1) - (|66I) leads 
to a set of equations in which the zeros determine the 
dynamical FPs. The only stable FP is found for wjj = 
/* finite. The corresponding values are 



w 



and w'j_ 

presented in Tab. [TTJ As a consequence we have v* = 0. 
We want to note that the static FP values of the two loop 
calculation in paper I and III have been used. We were 
interested in the non-asymptotic properties described by 
flow equations and since no real FP in statics is reached 
in two loop order we had to resum the static /3-functions 
in order to get real FP values @ . To each type of static 
FP (biconical or Heisenberg) two equivalent static fixed 
points exist differing in the signs of j± and 711 [TTj | . Ac- 
cordingly four equivalent dynamic fixed points exist with 
different signs in w'j_ and f±. They correspond to the 
directions of the external fields of the parallel and per- 
pendicular OP. 

The finite value of w\_ implies the relations 



£ + (r* x +Cx-2C = 0, (93) 



which follow from and (|66l) . The vanishing wjj leads 
to (p =0 as can immediately be seen from (IMI) . Using 
the first relation in and the third one, we obtain 



(94) 



Inserting the FP value of the static (-function (^ (see 
relation (105) in paper III) 



(95) 



(96) 



The dynamic critical exponents (|54[) in one loop order 
are therefore completely expressed in terms of the static 
exponents: 



c = 

v 2 

into the above equation one has 



Zj_ = Z ( J >± 



z <t>\\ 



(97) 
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These static exponents might also be taken from static 
experiments. All our numerical calculations are per- 
formed in d = 3 (e = 1). The numerical values of the 
static exponents <fi and v have been calculated in two loop 
order in paper I and are given there in Tab. Ill (v = v + 
therein) in two loop order resummed. In one loop or- 
der the two OPs have different dynamic critical expo- 
nents. Scaling is fulfilled only between the perpendicular 
OP and the secondary density. The parallel OP behaves 
like the van Hove model. This is demonstrated in FigUJ 
where the effective exponents defined in (f55|) - ([57)) have 
been calculated by using the flow equations in one loop 
order. At a flow parameter about I ~ e~ 15 for both, the 
biconical FP (solid lines) and the Heisenberg FP (dashed 
lines), the asymptotic values of the dynamic exponents 
z± and z m are reached. The classical value Z|| = 2, valid 
for both static fixed points, also is indicated by a straight 
line. The corresponding flow is presented in FigJ21 which 
proofs that the dynamic exponents in Fig|l]have reached 
their asymptotic behavior because the dynamic parame- 
ters are at their FP values at ln^ = —25. 



VIII. LIMITING BEHAVIOR OF THE 
DYNAMICAL (-FUNCTIONS IN 2-LOOP ORDER 
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The appearance of In v-terms in the two loop contribu- 
tion to the Cr^ -function, Eq ([70]) . changes the discussion 
of the fixed points considerably compared to the one loop 
case. In order to determine the dynamical fixed points 
of the current model in two loop order it is necessary 
to know something about the limiting behavior of the Q- 
f unctions. For this reason we will present the Q- functions 
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z ll 
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FIG. 1: Effective dynamic exponents at d = 3 calculated 
in one loop order using the one loop expression for the flow 
equations (|64|l , (|66|l . The effective exponents are calculated at 
the biconical FP (full lines) and at the Heisenberg FP (dashed 
line). Zii is valid for both FPs. 



FIG. 2: Flow of the parameters w' ± , s = w'^_/w'j_ and f± in 
one loop order at d = 3. The calculation has been performed 
for the biconical (solid lines) and the Heisenberg (dashed 
lines) FP. 



in cases where one or several dynamical parameters go to 
zero or infinity under definite conditions. This is neces- 
sary because some ^-functions exhibit singular behavior 
under these conditions, which influences the discussion 
of possible fixed points. It is anticipated that the critical 
exponents defined by the values of the C-functions at the 
FP are finite and real. 

fx) fx) 

The auxiliary functions X±, Lj_ , Ti, Ti and L\ , 

which appear in the ^-functions (|70|) and (|80|) behave sin- 
gularly in several limits of the parameters. Thus several 
FP values of the different parameters can be excluded 
due to diverging ^-functions. For a summary of the sub- 
sequent analysis of the £-function on the time scale ratios 
see Tab. UHl 

i) At first we will consider the two functions X± and 
L^ x ' in (|7i|) and (|7§j) . which appear in (r ± and depend on 
v only. These two functions remain regular if v grows to 
infinity. In this case one simply has X±(v — > oo) = 1 and 

L±(v — > oo) = 0. But for vanishing v both functions 
evolve a term proportional to In v. One gets 

Xj_(v -> 0) = ln(2w) , ii X) (w 0) = 1 - ln(2«) . 

Thus divergent In v terms appear in (j^ ( v 0) indepen- 
dent from the individual behavior of w± and w\\ because 
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FP 


u\ 


u\ 


< 


7y 


7-L* 


< 


W T 


/* 


B 


1.28745 


1.12769 


0.30129 


0.54201 


-0.17806 


1.55489 


=F 0.41958 


± 1.08563 


B 


1.28745 


1.12769 


0.30129 


-0.54201 


0.17806 


1.55489 


± 0.41958 


± 1.08563 


H 


1.00156 


1.00156 


1.00156 


0.85179 


-0.42590 


1.58136 


=F 1.38256 


± 1.24264 


H 


1.00156 


1.00156 


1.00156 


-0.85179 


0.42590 


1.58136 


± 1.38256 


± 1.24264 



TABLE II: FP values of couplings and timescale ratios for nu = 1, n± — 2 at d — 3. B indicates the biconical, H the Heisenberg 
FP. There are always two equivalent static FPs depending on the signs of the couplings 7. The FP values of the static couplings 
{m} and {7} are taken from the resummed two-loop results @], whereas w*_ and /*_ are calculated from the one-loop /3-functions. 
to! = v* = is valid in all cases. Corresponding to the two equivalent cases in statics and the sign of w"* there are equivalent 
dynamic FPs with corresponding signs of the FP value of the mode coupling /* . 



only the ratio v enters the function. 

ii) The dynamic ^-function (|80p of the parallel compo- 

(x) 

nent contains the three functions Ti, T2 and L,, defined 
in ([5T]). (|gg| and (JHU) which contain the ratio v. These 
functions, and therefore also £r» , remain non-divergent 
for vanishing v. One obtains 



Ti(v -> 0) = D ± I l+lnl 1 + — ) - — ln(l+v A 

v±J v±_ 

T 3 (v -*0) = -w* x D± I In I 1 + — J +ln(l + «x)) , 



L\*\v -> 0)=0. 
But they diverge when v is growing to infinity: 



(99) 



T\{v — > 00) = D± In 



T 2 (v -> 00) = -w* ± Dj_ In 
l\, x) (v -> 00) = 1 +ln 



1 

vx 



hi) 



(100) 



In contrast to the case i) the function in the paral- 
lel subspace evolves logarithmic terms lnz; in the limit 
v — > 00 and stays finite in the limit v — > 0. The above 
discussion is also independent of the individual behavior 
of w±_ and w\i because only the ratio v± stays always fi- 

nite and the three functions T±, T2 and L\, remain finite 
even for diverging time scale ratios if their prefactors are 
taken into account. 



Limit 



v - 


-+ 


~ \nv 


regular 


unaffected 


v — 


> 00 


regular 


~ Inv 


unaffected 


11)1 


-> oc 


~ lnwi 


regular 


12 


tW|| - 


-> 00 


regular 


~ lnw|i 


unaffected 


W X 


->■ 


regular 


~ lnu 


regular 


w \\ 





~ lnw 


regular 


unaffected 



TABLE III: Limiting behavior of the dynamic (^"-functions 



iii) Additional logarithmic singularities may arise in 
the dynamic ^-functions if the time scale ratios w± and 
w\\ grow individually to infinity independent of the be- 
havior of v. A closer examination of (1701) reveals that in 
the limit w± — > 00 the ^-function is proportional to 



Cr ± [w±_ -+ 00) ~ -7 ± In — 



(101) 



independent of the behavior of v. Quite analogously the 
same happens in (1801) when wu grows to infinity. One 
obtains 



Cr y (w\\ ~> 00) 



7 fln 



w \\ 

T 



(102) 



Supposing a finite (different from zero or infinity) FP 
value f]_ for the mode coupling parameter we may con- 
clude the following concerning the allowed FP values of 
the remaining parameters: 

a) From i) and ii) follows that v* has to be also different 
from zero or infinity, otherwise lnw contributions would 
lead to divergent ^-functions. 

b) From iii) follows that the finite v* only can be real- 
ized either by w\\ and w± both finite, or w\\ and w± both 
going to zero in the same way. The possibility that both 
time scale ratios are going to infinity in the same way is 
excluded because of the lnui|| and lnui^ terms appearing 
in this case. 



IX. GENERAL ASYMPTOTIC RELATIONS 

The FP values {en*} of the model parameters are found 
from the zeros of the /3-functions in Eqs. ([54]) - (I55|) . From 
the right hand side of the equations one obtains 



GO = 

, Cr,, ~ Ca ) = 



^7r Cr ^ 



= , 


(103) 


= , 


(104) 


= . 


(105) 



A FP which fulfills Eas.(fT03 |) -<fT05 j) has to be also a so- 
lution of 



= 



(106) 
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which follows from (|67[) . The £- function (r ± for the per- 
pendicular component of the OP relaxation is a complex 
function. Separating real and imaginary part leads to 



(107) 



As a consequence also the equations (|103|) for and 
(|106p for v are complex expressions. The £- function for 

r, is 



(108) 



with s defined in (|12ip . 

We anticipate that in a real physical system definite 
dynamical exponents exist, and therefore the dynamic £- 
functions have to be finite at the stable FP. As already 
mentioned in subsection lVIIIl the £- functions contain In v 
terms requiring a finite FP value v* in order to obtain 
finite dynamical exponents. Separating (|106[) into real 
and imaginary part one has 



l-k { /-l-k \ ll-k J- 

v (Cr , - Cr„ ) - v C: 



//*/-//* r\ 

v± — u > 
^ (r ± +v (C r± - CrJ = 0. 



From these two equations the FP relations 



Cf 



Cf' 







(109) 
(110) 



(111) 



immediately follow. The second equation in (|111[) implies 
that the dynamical exponent z$ ± in (|54|) can be written 
as 



z <f>± = 2 + Cf*L • 
From the first relation in (lllip follows 

= ■ 



(112) 



(113) 



This means that in the case of a finite FP value v* scaling 
between the OPs is valid. In order to obtain a critical 
behavior different from model C, the FP value /* of the 
mode coupling parameter also has to be finite and differ- 
ent from zero. Then from Eq. (|105[) follows 



(114) 



where the second relation of (II 1 1 p already has been used. 
Inserting ([55]) and ([931) into (II 14[) one obtains the relation 



Z <p± 1 Zl : 



(115) 



Dynamical strong scaling FP: 

In the case that w± and wu are finite at the FP, 



from ([103]) and (fT04[) the relation 



(116) 



is obtained, where (jlll[) already has been used. 
From ([54^) it follows immediately that the dynami- 
cal exponents have to fulfill the relations 



Hi 



(117) 



Thus in the case of strong scaling one dynamical 
exponent z exists only. The exact value of this ex- 
ponent can be found by inserting (j!17[) into (|115l) . 
One obtains 



v 



(118) 



(ii) Dynamical weak scaling FP: 

In the case that w± and tup are zero with v finite 
at the FP, Eqs.CEUS]) and (fTMj) arc trivially fulfilled 
and no additional relation between the ^-functions, 
and dynamical exponents respectively, arises. As a 
consequence two dynamical exponents exist. The 
first one 



Z <P» — z <p ± — z OP 



(119) 



for the OPs, follows from relation (|1 13|) . The sec- 
ond one 



z m — 2 zop 

v 



(120) 



for the secondary density, is obtained from (|115p . 



A closer examination of the /?- functions (|103l) - (|105[) . 
also with numerical methods in d = 3, reveals that no 
FP solution can be found where both, w± and tuii , are 
finite. Thus the only solution in d = 3 which remains is 
w'* = w'j* = iu| = with v* and w* finite. This result of 
course depends also on the specific numerical values [22[ 
of the static FPs (given in Tab [TTJ) used in the dynamical 
equations. The stable FP lies then in the subspace where 
the time scale ratios w\\ and w± approach zero in such 
a way that their ratios v and v± remain finite and in 
general complex quantities. In order to obtain the finite 
FP values for v and v± the two loop ^-functions may be 
reduced by setting w± and w\\ equal to zero by keeping 
their ratios finite. This will be performed in the following 
section. 



between the exponents. In summary, the condition that 
both v* and /* have to be finite leads to the two rela- 
tions (|113p and (I115[) between the exponents. Further 
relations are dependent whether the FP values of the 
time scale ratios w± and toy are finite or zero and lead 
to the following cases. 



X. CRITICAL BEHAVIOR IN THE 
ASYMPTOTIC SUBSPACE 

Since the asymmetric couplings j a always appear to- 
gether with the time scale ratios w a all terms propor- 
tional to these couplings drop out in the asymptotic limit 



13 



where w a — >• 0. It is convenient to introduce the real ra- 
tios 



g 



OH = El 

w' T' 



(121) 



Thus only s, q and f± remain as independent dynamical 
variables. 

The ratio s determines the behavior of the imaginary 
part of wj_ with respect to the real part, while the ratio q 
indicates the behavior of w\\ with respect to the real part 
of w j_. The complex parameters v± and v, introduced in 
(frJTj) and (|60|) , are expressed by s and q as 



1 



1 — is 



(122) 



in the following expressions. 



A. (-functions 

We discuss the behavior of the (-functions in the limit 
w± — > and w\\ — > for s and q constant. 
Case s ^ : 

For w± = and w\\ = the (- function (|70l) . reduces to 
C^ s) ({«}, *, g, /j.) = -^-{ 1 + |« ± (l («) 
+*_(*)*! (»)Li(«)) - i T ^-(2 a ;-Gs)L 1 (,s) 

-2L (s) + | + 4f ({«}, s, g) . (123) 

The functions cc-(s), xi(s), Lo(s), ^i( s ) an d Lr(s) are 
the same as in (JT5J) - dZTJ) with v± replaced by l[122jt . The 
same is true for (p ({ii},s,g), which has been defined 
in ([78| , and where also (|122[) has been used to replace v± 
and v. 

Performing the limit in the dynamical (-function (|80[) 
it reduces to 

(^({u},s, q )=(W({u},s, q ) (124) 

where ( r ^ ({it}, s, g) is the model A function (|83]) with 
relation (|122l) inserted into (|84|. 

Finally the function X2 in simplifies for vanishing 
time scale ratios to 



fx 



xr , (fi-) = f- 



(125) 



Inserting this expression into (|88[) and (|87j) . the dynam- 
ical (-function (|85l) reads 



Cf ) ({7} ) ^)=7i + i7f-^(l+^). (126) 



1 f 2 / f 2 



2 2 



The value of v± at the FP depends on how w' ± goes to 
zero in the critical limit I — > compared to u>". There 
are three possible scenarios: 

i) w" goes to zero faster than w' ± so that s — ► 0. Then 
u_L is turning to the real value 1. 

ii) and w'j_ behave in the same way so that the ratio 
s = so is constant. v± is in this case a complex constant 



1 + is 
1 - is 



(127) 



iii) goes to zero faster than w'[_ so that s — > 00. 
Then v± is turning to the real value —1. 

The third of the three scenarios above can be excluded 
from the discussion because some of the (-functions do 
not stay finite for v± = — 1. Finite (-functions at the 
FP and therefore well defined critical exponents may be 
obtained only in the first two scenarios. 

The (-function for scenario ii) are already given in 
(HUD - (gUD when (TT2Tj) is inserted. 

Case s — : 

For vj_ = 1 (s = 0) the (-functions (IT23)) and (fT24|) 
simplify to 



C^ 0) ({u},q,fx) = -fi 



l-^lnl-3 



+C£ UJ (M>«) ,(128) 

C ) (W>9)=C ) (M,9). (129) 
The model A functions (|78j) and (1531 are now 



CS 0) (W^) = ^(31ni-± 



4 1 



3 2 



(130) 



and 



C(W.«)-^(-i-5 



(131) 



where in L- X '(g), introduced in ([79]) and (|84|) . the rela- 
tions (|122[) with s = have been inserted. 



B. Fixed points in the asymptotic subspace 

Inserting the (-functions of the previous subsection 
into (|103l) - (|105l) one obtains the FP values in the asymp- 
totic subspace for s* finite, or s* =0. The results are 
presented in Tab II VI for the biconical (B) and the Heisen- 
berg (H) FP. It turns out that especially at the biconical 
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FP the values of the ratio q are extremely small, but def- 
initely not zero. Thus the asymptotic critical behavior 
in two loop order changes considerably compared to one 
loop (see section rvTTj) . Weak scaling as discussed in sec- 
tion |IX] is valid. The two order parameters scale with the 
same dynamic exponent zop from relation (|119|) . while 
the secondary density scales with a different dynamic ex- 
ponent z m given in (|120[) . The numerical values of these 
two dynamic exponents are also given in Tab II VI in two 
loop order. Note that the values of the dynamical expo- 
nents to the accuracy shown are independent wether the 
FP value of s is zero or not. 
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ZOP 
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2.18 
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F [10] 


0.83 









~ 1.5 


~ 1.5 


B 


1.232 


1.167 


1Q -86 
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1.131 


H 


1.211 


3.324 


• 10~ 8 





2.003 


1.542 


6 


1.232 


2.51 • 


10 -782 


0.705 


2.048 


1.131 


U 


1.211 


3.16 • 


1Q -66 


0.698 


2.003 


1.542 



TABLE IV: FP values of the mode coupling f± and the ratios 
q = w\\/w'x and s = w±/w'x in the subspace w\\ = 0, w± — 
and finite v — q/(l + is) for different cases of the biconical B 
and Heisenberg H FP in d — 3. For comparison results for 
model C and model F FPs are shown at n = 1 and n = 2, 
correspondingly. 



The comparison with the dynamical critical exponents 
in the cases when the OPs decouple statically and dy- 
namically into model C and model F shows the changes 
in the multicritical case where the exponents are changed 
but each component reflects the decoupled values accord- 
ingly. 
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FIG. 3: Effective dynamic exponents in the subspace w\\ = 
w± — with q and s finite in d — 3. The static values 
are taken for the Heisenberg FP and for the biconical FP. 
The non- asymptotic region is extended by a factor 10 at the 
biconical FP. For the static FP values see Tab. [Ill f° r the 
dynamic FP values Tab. IIVI 



C. Effective exponents in the asymptotic subspace 

The flow of the parameters q, s and f± can be found 
by solving the equations 



ds 

a 

,df± f±f Aas) 



2Cm + &[C£ a) ] 



(132) 
(133) 

(134) 



The (^-functions in the above flow equations are the re- 
duced expressions (|123l) . (|124[) and (|126|) , which are func- 
tions oiq,s,f±. We consider the case s ^ since the FP 
s* = is reached only starting with s = 0. From the 
solution of Eqs. flll^ -flTM )) the flow q(l), s(l), fj_(l) is 
obtained, which is used to calculate asymptotic effective 



dynamic exponents 

*i">(0 = 2 + »[c& ) (?(i),«(0,/x(0)] 

-sm[4^(q(l) lS (l)J ± (l))} , (135) 

(136) 
(137) 



~! as) (0 = 2 + C^ ) (q(l),s(l)) , 



(as) . 



t s) (l) = 2 + Cr(q(l),s(l)J ± (l)) 



They can be calculated for different static fixed points, 
i.e. biconical or Heisenberg FP, as presented in Fig|3] 
The values of well as 7^_, 7|, used in 

the current calculations can be found in Tab. II. At both 
fixed points weak scaling, as discussed in section IIXI is 
fulfilled. The difference to the one loop result is now 
that the dynamic exponents z± and zn of the OPs are 
equal in the asymptotic region, while z m stays different. 
Moreover the transient exponents in two loop order are 
very small compared to one loop. There the effective ex- 
ponents reach their asymptotic values about I ~ e~ 15 as 
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can be seen from FigJTJ In two loop order the asymptotic 
region is of magnitudes smaller. From Fig(3]one can see 
that at the Heisenberg FP the flow parameter has to be 
of the order / ~ e~ 1500 to obtain the asymptotic values 
of the dynamic exponents. At the biconical FP I has to 
be even of the order I ~ e _1500 ° (note that there is a fac- 
tor 10 between the cc-scales in Fig|3J| to reach asymptotic 
values. However as will be seen in the next section the 
subspace will not be reached by the flow in the complete 
parameter space for reasonable values of I. 



XL GENERAL FLOW AND 
PSEUDO-ASYMPTOTICS 

Although in general the dynamic flow equations have 
to be solved in the full parameter space, the results for the 
effective exponents presented in the previous subsection 
are obtained from the flow equations which already have 
been reduced to the subspace w± = w\\ =0. The reason 
to do this is that the flow and the ^-functions in the full 
parameter space shows some peculiar behavior. 

In the non-asymptotic region the flow is generated by 
the system of equations for four parameters which are 
w' ± , s, w\\ and f± obtained from Eqs. (|64l) - (|66|) and 
([1211) . The static parameters are taken at their FP val- 
ues given in Tab. |TTJ Due to the presence of the static 
asymmetric couplings ji and the mode coupling f± an 
imaginary part of tui is produced even if one starts with 
a zero initial value. Starting with a typical set of initial 
values, i.e. w\\Qo) = 0.3, w' ± (lo) — 0.6, s(Iq) = 0.5 and 
f±(lo) = 0.4 at the flow parameter value ln/o = — 1, the 
effective exponents in the complete parameter space have 
been calculated in d = 3. The result is presented in FigfJ] 
for both static fixed points, where the solid lines are the 
results of the two loop calculation and the dashed line 
is a result of a complete (flow and effective exponent) 
one loop calculation. However the static FP values from 
Tab. [II] have been used also in the dynamic one loop 
flow. There it seems that in two loop order the same 
results as in the one loop calculation are obtained. z± 
and zffi are getting close together (solid lines) for flow 
parameters I < e~ 100 and seem to coincide even numer- 
ically with the corresponding results in one loop order 
(dashed line). This is the type of weak scaling in one 
loop order, which also can be seen from Figfl] and in 
qualitative contradiction to the discussion in the previous 
sections (see FigOJ) - But the examination of the flow of 
the dynamic parameters reveals a fundamental difference 
between the one and two loop calculation. In one loop 
order the dynamic parameters w' ± , s, Wu and f± merge 
to the FP values when the effective exponents turns over 
in their constant asymptotic values, which is presented 
in Fig[2] This happens in the region about / < e~ 15 and 
the dynamic parameters stay constant for all lower flow 
parameter values. In two loop order the situation is dif- 
ferent. Although the two loop results for the effective 
exponents look like one has reached the asymptotic re- 
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FIG. 4: Effective dynamic exponents in the background us- 
ing the flow equations (|64l) . (|66[) in two loop order in d = 3 
in the complete dynamical parameter space (full lines). For 
comparison the effective dynamic exponent 2^' in one loop 
order is shown (dashed line). 



gion (the exponents seem to be constant), the dynamic 
parameters in contrast are far from their asymptotic FP 
values. This is presented in Fig (5] The parameters w'j_ 
and s are still increasing and obviously have not reached 
a FP value. At the first glance the flow of f± seems to 
have reached a FP value (see lowest plot in FigH]). But 
a closer examination shows that this is not the case. f± 
is constantly increasing with a very small slope as can be 
seen from the inserted small figure, where both axes has 
been enlarged. Actually the set of two loop /3-functions 
does not have a zero for finite w± and Wu , and therefore 
no FP exists in the parameter region of FigJS] Thus the 
effective exponents in two loop order in Figf?] only show 
a pseudo-asymptotic behavior completely different from 
real asymptotic behavior (there z± and z?, have to 
be equal) discussed in section IXCl (see Fig. [3]). Even if 
one draws the x-axis in FigfJ] down to In / = —20000, as 
done for the flow in FigfSJ the picture remains to be the 
same, that is the effective exponents seem to be constant, 
with the exception that the background behavior is not 
longer visible because the region is too small. Also if one 
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FIG. 5: Flow of the parameters w' ± , s and f± in the full 
parameter space 



changes initial conditions of the parameters the qualita- 
tive result remains the same. The different flows merge 
within a region of InZ = —150 to the same result. 

In order to get some insight how this pseudo- 
asymptotic behavior is possible, in Fig|5] the relative 
slopes 



1 don 
on d In I 



a; 



(138) 



for the parameters on chosen to be w'j_ and f± have been 
calculated. One can see that the relative changes in these 
parameters drop down to very small values. As a first 
consequence one has to calculate down to extremely small 
flow parameter values InZ < — 10 8 where one can expect 
to leave the pseudo-asymptotic region. But although the 
/3-functions cannot be zero in the considered parameter 
range they may reach values which are so small that they 
cannot longer be separated numerically from zero. This 
means that coming from the background it is impossible 
to pass the pseudo-asymptotic region numerically into 
the real asymptotic region. This is the reason why in the 
previous section the flow has been started in an asymp- 
totic subspace. 

Due to the presence of logarithmic terms in the 
timescale ratio v the FP value of v has to be finite. As 
follows from Tab. IIVI via Eq. (I122[) it turns out to be 



very small leading to very large (negative) values of In v. 
So one expects that the In u-terms begin to dominate Cr 
in a certain region of In I near the asymptotics. Making 
the In i>-terms explicit one may rewrite (j^ , given in Eq. 
([?D1l. as 



1 7||Z?j 



1 l\\Dj 



u 2 

"x 1 

m v 

36 



2 1 + wa V 3 2 1 
remaining terms . 



wj_ 



In v 



(139) 



Inserting Eq. (|122[) . the essential term is the real part V 
of the prefactor of In q. One obtains 



Cr^ =V']nq 
The prefactor V' is 
2 1 



remaining terms . 



v - 

with 



36 



IP 2 



A 



w' ± [(l + w' ± + w' ± s 2 )j± - sF] 
(1 + w' ± ) 2 + {w'j_s) 2 



B = 



w' ± sj ± - (I + w' ± )F] 
(1 +w' A ) 2 + (w'^s) 2 



(140) 



(141) 



(142) 



(143) 



In FigJTJthe behavior of the In g-contributions to (p ± at 
the biconical FP is presented. Although the In q term al- 
ready reaches very large negative values (q is very small) , 
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FIG. 7: Contribution of the V'lnq term in (fT40)l (lower fig- 
ure). It is of the same magnitude as the other terms in the 
(^-function. V' and In q have the same behavior with opposite 
exponents. This is demonstrated in the upper figure, where 
the decadic logarithm of V' (solid curve) and the negative 
decadic logarithm of — \nq has been drawn (dashed curves). 

as expected, this is compensated by the prefactor V 
which has very small values in the considered region. In 
the upper part of Fig|7] the decadic logarithm of V and 
— In q have been plotted. Both curves show a similar be- 
havior and a small difference. In the lower part V In q is 
calculated from (|141[) - (|143[) . As a consequence of the 
results in the upper part of the figure the numerical val- 
ues are about —0.5 and one can see that the term is far- 
away from being the leading one. There are other con- 
tributions to Cr_L which have the same magnitude. Thus 
the ^-function is not in the asymptotic region as has been 
also indicated by the flow in FigJSJ 

Thus one expects in the experimentally accessible re- 
gion non-universal effective dynamical critical behavior. 
This is described in the crossover region to the back- 
ground by the flow equations together with a suitable 
matching condition related to the temperature distance, 
the wave vector modulus etc. The initial conditions have 
to be found by comparison with experiment. 



XII. CONCLUSION 

Our two loop calculation for the dynamics at the mul- 
ticritical point in anisotropic antiferromagnets in an ex- 
ternal magnetic field leads to a FP where the OPs charac- 



terizing the parallel and perpendicular ordering with re- 
spect to the external field scale in the same way (strong 
dynamic scaling). This holds independent wether the 
Heisenberg FP or the biconical FP in statics is the stable 
one. The non-asymptotic analysis of the dynamic flow 
equations show that due to cancelation effects the criti- 
cal behavior is described - in distances from the critical 
point accessible to experiments - by the critical behav- 
ior qualitatively found in one loop order. That means 
the time scales of the two OP components become al- 
most constant in a so called pseudo-asymptotic region 
and scale differently. 

So far we have not included the non-asymptotic flow 
of the static parameters which are expected to lead to 
minor deviations from the overall picture. Another item 
would be the study of the decoupled FP since in the non- 
asymptotic region the OPs remain statically and dynam- 
ically coupled and the behavior depend on the stability 
exponents how fast these effects decay. This in turn de- 
pends on the distance of the system in dimensional space 
and the space of the OP components from the stability 
border line to other FPs than the decoupled FP (see Fig. 
1 in paper I). 

The numerical results presented in this series of papers 
have been calculated for dimension d = 3. One might 
speculate that the peculiar behavior found is specific to 
the dimension of the physical space rather than to the 
multicritical character of the specific point. This aspect 
was out of the scope of this series of papers. We note that 
the two critical lines meet at the multicritical point (bi- 
critical or tetracritical) tangential. This has been taken 
into account for the nonsaymptotic behavior by choosing 
a path approaching the multicritical point without meet- 
ing one of the two critical lines [23[. The nonasymptotic 
behavior in fact is more complicated since two critical 
length scales are present in the system. This has to be 
taken into account when studying the crossover behavior 
in approaching one of the critical lines [24| . 

Only recently a bicritical point has been identified by 
computer simulation [2o| . The corresponding FP has 
been identified as the Heisenberg FP which corresponds 
to the type of phase diagram obtained. It seems to be 
difficult to look for situations where a phase diagram 
containing a tetracritical point is present. Even more 
complicated would it be to identify the dynamical char- 
acteristic of this multicritical point, where - coming from 
the disordered phase - two lines belonging to different 
dynamic universality classes meet. The dynamical uni- 
versality class of the case with a n = 1 OP (model C) has 
been studied in [26|, [27] with different results leading to 
critical exponents larger than expected. The dynamical 
universality class of the case with a n — 2 OP (model F) 
case has been studied by computer simulations in [28l [29| . 
The methods of these simulations might be extended in 
order to be used also in the case of the multicritical point 
studied in this paper. 
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Appendix A: Calculation of the dynamic vertex 
functions of the OPs 

In perturbation expansion up to two loop order the 
functions and fi^ ||5 which appear in 

and (|28|) . can be written as 

+«S1&->*II>*> W )> (Al) 



The dynamic propagators a j_ and a!« are defined as 

a' x = f ±(ei 2 + (fc + fc') 2 ) + M' 2 , (A9) 
a\ = f || {q 2 + {k+k'f) + \k' 2 . (A10) 
The two loop contributions to the dynamic vertex func- 
tion of the orthogonal components (IA1[) and (| A2|) have 
the structure 

ftS^(£L,C B ,*,«) = §f J .^W'tJi>(e 1 . I *,a,) 

+ (f±7±-iff)7J-^+(a,e||,*,w) (AH) 



(A2) 



(A3) 



The superscript (iL) indicates the loop order. Of course 
all functions considered depend on all model parameters 
(couplings and kinetic coefficients), but only the indepen- 
dent lengths £xj £|b k and ui will be mentioned explicitly 
in the following. The one loop contributions are 



(l^l (U,k, w) = (f±7x - iff)7±J±(£±, w) » (A4) 



and 



0«.- ■/•■-•) = -^rxfixisJ%;^(a,fe,w) 

r x« x is^+ (£x,£||,fc,w) 



2 , 



(T3j 



6 



rx7± - (Cx,C|| , fc, w) -(A12) 



Note that both two loop functions differ only in terms 
containing the static fourth order couplings u± and u x - 
The remaining contributions are the same in both func- 
tions apart from a factor j± and ig respectively. The 
function F^ + is defined as 



Gj|i(a,*,w) = (fx7x-i5)i^x(C||,fc,w) > (A5) F^ + (U,Z h k,u>) = (fx7x - igfP^faM 



rt S?,(fD.*.«) = Wll«ll.*' w ) 
The one loop integrals Jx and Jm in (|A4[) - (|A6[) read 

1 



(A6) 



Jk 

I\\(€\\,k,u 



7T> (A?) 



fc' (G 2 + (fc + fc') 2 )(-i^ + al) 

/ -o —s • (A8) 

k> (C 2 + (k + k>r)(-iu; + a>) 



-I ( F {T ** } (a , £|| , - 7|| F^?* ] (U , C|| , k, u) 



+F 



The first two loop contributions in (IA1 1[) come from the 
bicritical model A. The integrals W , 7 '} are defined by 



19 



W$i\t±,k,u>) [ I ( G 2 +(A+jfe0 a) (G 2 + ^ ) ( fI 2 +(feW)2)( _ ki;+A _ L _ L+ _ L) . ( A14 ) 



k' Jk" 



(e: 2 +(fc+fco 2 ) +fe" 2 ) (q 2 +(k'+k") 2 ) (-iw+^xu ii ) 



i 



with The further two loop contributions in (|Allj) - (|Al3j) are 

marked with superscripts (Ti), which indicate the diffcr- 
A±± + ± = f ±(ei 2 +(fc+fc') 2 ) + (Cl 2 +fc" 2 ) ent graph topologies. The explicit expressions are 

+f ± (ei 2 + (fc'+fc") 2 ) (A16) 



and 



I 

+f||(^ 2 +(fc'+fc") 2 ). (A17) 



^ + v ; (ei 2 +(fc+/c') 2 )(-^+al)(-^+^±± + ±) VC + fc " 2 £i 2 +(fe'+fc") 2 y 

F {T -iL x) (£x,£\\,k,u))= [ ( 7—, r „ 2r " 7 " , (A19) 

^ + k " ' Jk>Jk» (G 2 + (fc+fc') 2 )(^ 2 +fc" 2 )(-i^+«j|)(-i^+^||||) ^ ^ 

Ff^f^fc.w) = / / — -i , (A20) 

^ + V ' i fc J fc »(ei 2 + (fc + fc' + fc") 2 )(-i^ + a ± ) 2 (-iw + /3i) V ' 



p{T^) (f ,s[ f 7± Afc' 2 -i g [(fc'+fc") 2 -fc" 2 ] T ± j ± -ig r+7±+ig \ fA21 s 

^ + lU ' ' j AJ fc » (C 2 +(fc+fc') 2 )(-iu;+al)2(- 1 c+A ±±+± )Ul 2 +fc" 2 ei 2 +(fc'+fc") 2 / ' 



. m , /• /• 2rn7 2 Afc' 2 

F ( -l*'(£±,Z\\,k,e>) = / / ^5 s 5 — , (A22) 



£(T6x),, fc ^ = /" y 7± Afc" 2 +i g [(fc+^+fc") 2 -(fc+fc0 2 )] / r ± 7±-jg r+ 7± +i g 

^ + lU ' ' WJ Jk' Jk4^ 2 Hk+k') 2 )(^+a' ± )(-ioj+a'[)(-iu J +S ±±+± ) \!^ 2 +(k+k>+k") 2 ^ 2 + (k+k")\ 

f±j±-ig ( r ± 7j_-ig 7j_ 



fe < 7 fc »(ei 2 +(fc+fc') 2 )(-i^+aT)(- iw +^) V - iLJ+a '± ei 2 + (fc+fc'+A:") 2 
7± AA:" 2 +i.g[(A:+fc'+fc") 2 -(fc+fc') 2 )]\ 



(A23) 



(£i 2 +(fc+fc'+fc") 2 )(-i^+a / ± ) y 

with the dynamic propagators 

± = f ± (ei 2 +(fc+fc'+fc") 2 ) + A(fc' 2 +fc" 2 ) , (A24) 
5 ±±+± = t ± (ei 2 +(fc+fc') 2 ) + f 1 (G 2 +(fc+fc'+fc") 2 ) + f jl (£I 2 +(fc+fc") 2 ) (A25) 
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which are both invariant under an interchange of k' and k" . 

The two loop contributions to the dynamic vertex function of the parallel component (|A3[) has the structure 



" |f||"xT||^|; ) (^,e||,fc,^) + f || ^^ ||<?|| (^,Cn,fc,^). (A26) 



The function F, i is defined as 

<t>\\n 



The integrals W-^' 4 -'' in the two loop contributions from the bicritical model A are 

011*11 

Mw)r - r r i 



w hh^ M JJ k „ (e ir 2 +(fc+fco 2 )(^ 2 +fc" 2 )(ef 2 +(^") 2 )(-i-+^ ) ' (A28) 

KTh^^^ = L L (^ 2 +(fc+fc') 2 )(G 2 +fc" 2 )(G 2 +(fc'+fc") 2 )(-^+^ii i+i ) (A29) 

with the propagators 

A = f || (^ 2 + (k+k') 2 ) + f || (^ 2 +fc" 2 ) + f || (£- 2 + (fc'+n 2 ) (A30) 

and 

A [[x+± = f||(^ 2 +(fe+fc') 2 ) +f+(ei 2 +fc" 2 ) +f ± (G 2 +(A : '+A : ") 2 ) . (A31) 
The remaining two loop contributions in (|A26j) and (|A27|) are 

F iT3 j\t h k,u)= [ [ _ , - 2f|1 ^" , (A32) 

" A' A» (^ 2 +(A:+fc') 2 )(^ 2 +A:" 2 )(-i W +aj | )(-ia;+A )' V ' 



f±7_L-i3 f+7 ± +i.g 



A' A» (^ 2 +(fc+A:') 2 )(-^+«||)(-^||± + ±) Ul 2 + fc " 2 + G 2 + (fc'+fc") 2 J ' (A33) 

i^ (T4||) (£l,fc, W ) = / / „ , (A34) 

' i fc J fe "(^ 2 + (/c+/c'+fc") 2 )(-i^ + aj | ) 2 (-i W + /3||) ^ j 



V " ' A' A" (^ 2 +(fc+fc') 2 )(^ 2 +fc" 2 )(-i^+«I|) 2 (-i^+^ ) 

m^M=f I 1 3 ^ Afc, V i ^( fc,+fc ") 2 - fc " 2 ] f^g + ), (A36) 



fe < A» (^ 2 +(fc+fc') 2 )(-^+«||) 2 (-^+^n± + i) VG 2 + fc " 2 G 2 +(fc'+fc") 2 

f||7 2 Afc" 2 / 1 1 



^U'"'"'' J k J k ,,(q 2 + (k+k') 2 ){-iu+a\ l )(-iLj+a'[)(-ioj+S m ) \^ 2 + (k+k'+k") 2 ' £~ 2 +(fc+fc") 2 

;,2 



r H^ f r ll , 1 /\ , 



l k Jk>iq 2 + (k+k') 2 )(-iLo+a' l l)(-iLj+p ll ) y-iuj+a'^ q 2 + (k+k>+k") 2 \ -iw+aj| 



(A37) 



21 



The additional dynamic propagators are 

S =f||(^ 2 + (fc+A:') 2 ) +f | |(^ 2 +(fc+fc'+fc") 2 ) +fn(^ 2 + (fc+fc") 2 ) 

and 

/3,l = f || (q 2 +(k+k'+k") 2 ) + X(k' 2 +k" 2 ) . (A38) 

The integrals contained in (|A14|) - (|A23|) and (|A28I) - (|A37|) are of the same type as already has been presented in 
[30I ] (see Eqs.(A19) - (A26) in the appendix therein). The e-poles of these integrals can be found in Eqs.(C2) - (C9) 
of the same reference. 



Appendix B: Dynamic Z-factors in two loop order 



Within the minimal subtraction scheme of the renormalization group calculation one has to collect in two loop order 
the pole terms of order 1/e 2 and 1/e in the functions and f in (|27[) . The resulting dynamic renormalization 



factors are 
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The coupling D± and the functions A±, and and Xj_ are defined in ([7T|) -([74 |l . The pole terms of the function 

T^f\ + are collected in the renormalization factor 
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The functions T x , T 2 and T A have been introduced in (|8TT) - (|84|) . 

The renormalization factor Z\ is identical to the one of model F [21] with all parameters of the perpendicular 
subsystem. With Q defined in (|88[) . one gets 



Zx ~ 1 ~ e 2 1 1 + 2 4e w' 
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